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THE  APPLICATION  OF  CERTAIN  PROJECTIONS  TO  THE 
CONSTRUCTION  OF  GENERAL  MAPS  OF  CANADA 


By  Louis  B.  Stewart,  Professor  of  Surveying  and  Geodesy 

For  the  representation  of  a  considerable  portion  of  the  surface  of  the 
earth  numerous  methods  have  been  devised  and  used  in  the  construction 
of  maps,  and  among  the  best  of  these  are  to  be  included  the  different 
modifications  of  the  conical  projection  with  two  standard  parallels.  It 
is  here  proposed  to  examine  some  of  these  modes  of  representation  with 
a  view  to  determining  their  adaptability  to  maps  of  the  Dominion  of 
Canada. 

In  all  the  modifications  of  the  projection  to  be  examined  the  parallels 
of  latitude  are  represented  by  arcs  of  circles  having  a  common  centre, 
and  the  meridians  by  straight  lines  radiating  from  the  same  centre. 
Also,  the  scale  is  true  along  the  two  standard  parallels,  situated  near  the 
northern  and  southern  limits  of  the  map.  The  scale  in  longitude  is 
therefore  slightly  reduced  between  the  standard  parallels,  and  enlarged 
outside  of  them.  Various  additional  assumptions  may  be  made,  giving 
maps  suited  for  special  purposes.  Thus,  it  may  be  assumed  that  areas 
shall  be  correctly  shown;  or,  that  all  small  portions  shall  be  shown  in 
their  true  forms. 

An  obvious  advantage  of  these  projections  is  that  they  may  be 
easily  constructed.  If  the  scale  of  the  map  is  not  too  large  the  parallels 
may  be  drawn  by  the  aid  of  beam  compasses;  and,  the  degrees  of  longi¬ 
tude  having  been  laid  off  along  the  standard  parallels,  the  meridians 
may  then  be  drawn  by  joining  the  points  of  division  by  straight  lines. 
If  the  scale  is  large  the  standard  parallels  may  be  constructed  by  means 
of  the  rectangular  coordinates  of  their  points  of  intersection  with  the 
meridians,  the  axes  of  coordinates  being  the  central  meridian  and  lines 
at  right  angles  to  it  through  the  points  of  intersection  of  the  parallels. 

Other  advantages  are:  the  meridians  and  parallels  intersect  at  right 
angles;  the  scale  is  the  same  at  all  points  on  a  parallel,  either  in  latitude 
or  longitude,  and  therefore  the  scale  error,  or  the  distortion,  does  not 
increase  with  increase  of  longitude  from  the  centre  of  the  map,  as  is  the 
case  with  some  projections  that  have  been  extensively  used.  These 
projections  are  therefore  well  adapted  for  maps  covering  a  wide  range 
in  longitude,  and  at  the  same  time  having  a  considerable  extent  in 
latitude. 
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For  the  sake  of  comparison  brief  descriptions  of  the  Polyconic  and 
Bonne’s  projections  will  be  given  at  the  end  of  this  account. 

The  following  modifications  of  the  type  projection  will  be  considered 
in  turn: 

(1)  The  scale  in  latitude  is  to  be  true  at  all  points,  and  the  standard 
parallels  are  to  be  selected  arbitrarily. 

(2)  The  scale  in  latitude  is  to  be  true  at  all  points,  and  the  maximum 
scale  error  in  longitude  between  the  standard  parallels  is  to  be  equal, 
with  opposite  sign,  to  that  of  either  extreme  parallel. 

(3)  The  area  of  any  portion  of  the  map  is  to  be  correctly  shown. 
This  projection  is  known  as  Albers’  Equivalent,  or  Equal-Area,  Pro¬ 
jection. 

(4)  All  small  portions  of  the  map  shall  be  shown  in  their  true  forms. 
This  gives  Lambert’s  Conformal  Projection. 

In  order  to  institute  a  comparison  between  these  projections  the 
following  quantities  will  be  determined  in  each  case: 

The  scale  errors  in  latitude  and  longitude; 

The  error  in  the  representation  of  areas; 

The  error  in  the  representation  of  directions  at  a  point. 

The  following  notation  will  be  used: 

0  will  denote  the  latitude  of  a  parallel; 

0i  and  02  those  of  the  southern  and  northern  standard  parallels; 

0'  and  0"  those  of  the  extreme  parallels; 

co  the  longitude  of  a  meridian;  reckoned  from  the  central  meridian 
of  the  map; 

6  the  angle  at  the  centre  of  a  developed  parallel  corresponding  to  a 
given  value  of  co; 

r  the  radius  of  a  developed  parallel; 

n  the  ratio  6  :  co,  known  as  the  “constant  of  the  cone’’; 

p  the  radius  of  curvature  of  the  meridian  at  a  point  in  given 
latitude; 

N  the  length  of  the  normal  at  a  point,  terminated  in  the  axis  of 
rotation  of  the  terrestrial  spheroid; 

P  the  radius  of  a  parallel  of  latitude; 

m  the  length  of  a  given  meridian  arc. 

Convenient  expressions  for  p  N  and  P  are  as  follows: 

a(l-C)  (1) 

(1  —  e2  sin2  0)f 


N  = 


a 

(1  —  e2  sin2  0)  J 


P  =  N  cos  0 


(3) 
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in  which  a  denotes  the  radius  of  the  equator,  and  e  the  eccentricity  of  the 
elliptic  meridian.  Substituting 


these  become 


sin  \[/  =  e  sin  <f> 

(4) 

P  =  a{  1  —  e2)  sec3  \f/ 

(5) 

N  =  a  sec  \J/ 

(6) 

P  =  a  sec  i p  cos  4> 

(7) 

For  the  Clarke  spheroid  of  1866: 

log  a  (in  miles)  =3.5980536 
log  e  =2.9152513 

log  (1-c2)  =1.9970504 

The  length  of  a  meridian  arc  may  be  found  by  means  of  the  expres¬ 
sion  : 

m  =  [1.839 19449]A</>  (in  degrees) 

—  [1.3043545]  cos  2 </>0  sin  A </> 

+  [2.3301509]  cos  4 <£0  sin  2A<£ 

—  [5.45093]  cos  60o  sin  3A</>  +  (8) 

in  which — 

A (f>  denotes  the  difference  of  latitude  of  the  extremities  of  the 
arc,  and 

</>o  the  mean  of  the  extreme  latitudes. 

The  numbers  in  rectangular  brackets  are  the  logarithms  of  numerical 
coefficients. 

For  a  short  meridian  arc — less  than  1° — the  following  expression  may 
be  used: 


m  =  p.A</>  sin  1 


n 


(9) 


A0  being  the  amplitude  of  the  arc  in  seconds,  and  p  the  radius  of  curva¬ 
ture  at  its  middle  point. 

The  Supplement  to  the  Manual  of  Instructions  for  the  Survey  of 
Dominion  Lands  contains  a  table  giving,  for  every  10'  of  latitude  Irom 
42°  to  70°,  the  values  of  the  logarithms  of  p  sin  1"  (p  being  there  denoted 
by  R)  N  sin  1"  and  P  sin  1",  the  unit  being  the  Gunter’s  chain.  The 
log.  of  either  of  these  quantities  in  miles  may  be  found  by  adding  to  the 
tabular  quantity 

log - - - =  3.4113351. 

80  sin  1" 

The  table  on  p.  11  below,  giving  the  lengths  of  degrees  of  latitude  from 
40°  to  70°,  may  be  used  for  finding  the  length  of  a  long  arc,  the  fractional 
degrees  being  found  by  eq.  (9). 

We  shall  proceed  to  the  investigation  of  the  projections  enumerated 
above. 
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(1)  The  scale  in  latitude  is  true  at  all  points,  and  the  standard 
parallels  are  selected  arbitrarily. 


F/g  / 

PffojEC  r/o/v  No  2  ' 

Scale  _  /  Inch  =  500  Miles 


In  Fig.  1  AA'  and  BB '  are  the  standard  parallels,  their  radii  being 
T\  and  f2,  P\  and  denoting  their  true  radii.  Then  we  have 

VA  _AA'  _  N\  cos  </>i .  co  Pi 
VB  BB'  N2  cos  (f> o.o)  P2 
VA  Pl 

VA-VB  Pi- P2 
!n  Pi 


or 


m  Pi  — Po 

PlW 


or 


r 


P1-P2 


(10) 
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m  denoting  the  length  of  the  meridian  arc  AB.  Also 

r2  =  Y\  —  m 

Similarly,  the  radius  of  any  developed  parallel  is 

r  =  rizhm' 


m’  being  the  length  of  the  corresponding  meridian  arc. 
Again, 

AA'  =  r±6  =  P{to 

so  that 


6  _P\  _P i  — Pi 

co  T\  m 


(11) 

(12) 


(13) 


Also,  denoting  the  scale  ratios  in  latitude  and  longitude  by  a  and 
respectively,  we  have 

(7=1 


, _  rO  _  r  _P\—  Pi  r 
Poj  P  m  P 

Also  the  area  scale  is  given  by  the  expression 

rdrdd 

pd4>  .  Pdo: 

But  dr  =  pd(})1  and  dd  =  ndu,  this  expression  becomes 


(14) 


so  that  <j'  is  also  the  area  scale,  which  is  otherwise  evident. 

To  investigate  the  local  misrepresentation  of  direction  at  a  point 
let  us  consider  an  infinitesimal  portion  of  a  line  drawn  upon  the  surface 
of  the  earth  and  having  the  azimuth  a,  and  assume  that  a  meridian  is 
drawn  through  one  extremity,  and  a  line  perpendicular  to  the  meridian 
through  the  other  extremity,  thus  forming  a  small  triangle,  the  lengths  of 
whose  sides  are  pd(f>  and  Pdoi,  respectively;  so  that  we  have 


tan  a 


Pdu ). 
pd4> 


Then,  if  a  denote  the  angle  on  the  map  corresponding  to  a,  we  have 

/ 

<T 


tan  a' = 


•'Pdoi  o 


=  —  tan  a 
apd(f)  a 


(15) 


An  expression  for  a'  — a  may  be  obtained  in  the  form  of  a  series  Eq. 
(15)  is  of  the  form 

tan  x  —  p  tan  y, 

which  by  expansion  in  series  gives 

x  —  y  =  q  sin  2 y  +  §52  sin  4 y  +  fg3  sin  6y 

in  which 

P-1 
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Applying  this,  we  have 


<7 


t  \  9 

<T  ~  <7 


a  —  a  = 


sin  2 a  +  4  (  - - -  )  sin  4a + 

cr'  +  cr  Vcr'  +  cr 


If  a  =  45°  the  first  term  of  this  series  has  its  maximum  value,  the  second 
term  vanishes,  and  the  third — as  well  as  each  succeeding  term  that 
does  not  vanish — is  extremely  small,  so  that  the  maximum  value  o 
a'  — a  in  seconds  is,  very  approximately: 


a  —  a 


a  —  a  = 


(a'  +  a)  sin  1 
In  the  present  case  this  becomes 

<x'-l 


tr 


a  —  a  = 


rt 


(16) 


(17) 


(a'+l)  sin  T 

In  order  to  determine  the  latitude,  between  those  of  the  standard 
parallels,  for  which  the  scale  error  in  longitude  is  a  maximum  we  have 


da'  n  d 

—  =  0,  or  n  — 
dp  d(f) 


-U)=°. 

i>\Pj 


But 


d_ 

d(f)  \  i'/  p 

the  condition  for  a  maximum  is: 


dr  d$ 

/  \  p— ~r — 

^ r  j  d$  dP^  ; 


Then,  as 


Also 


dr  dP 

P  —  =  r  — 

d<[)  d<t> 


r  =  ri  —  m 
dr  dm 


=  —  P- 

i 


d  <p  dp 

P  =  a  cos  4>  (1  —  e2  sin2 

dP  m  _3 

"  =ae2  sin  p  cos2  0(1  —  e2  sin2  0)  2 


dp 


—  a  sin  0  (1  —  e2  sin2  0)  2 
ae2  sin  0  cos2  p  —  a  sin  0(1  —  e2  sin2  p) 


(1  —  e2  sin2  0)  2 
a(l  —  e2)  sin  p 

(1  —  e2  sin2  p)% 

=  —  p  sin  0 

the  condition  for  a  maximum  becomes 

pN  cos  0  —  (r  1  m)p  sin  p 
N  N. 


or 


tan  0  = 


(IS) 


r\  —m 

This  equation  must  be  solved  by  trials. 
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The  coordinates  of  the  intersections  of  the  meridians  and  parallel 
may  be  computed  by  means  of  the  equations: 

d  =  no3  x  =  r  sin  6  (19) 

y  =  r(  1  —  cos  0) 

The  following  table  gives  the  length  in  miles  of  1°  of  latitude  for 
different  latitudes. 


Lat. 

Length  of  1° 

Lat. 

Length  of 

O 

m 

O 

m 

40 

68.99917 

55 

69 . 18039 

41 

69.01131 

56 

69.19181 

42 

69.02351 

57 

69.20306 

43 

69 . 03575 

58 

69.21414 

44 

69.04801 

59 

69.22502 

45 

69 . 06029 

60 

69.23570 

46 

69.07256 

61 

69.24616 

47 

69 . 08482 

62 

69.25639 

4S 

69 . 09704 

63 

69 . 26637 

49 

69.10921 

64 

69.27609 

50 

69.12132 

65 

69 . 28555 

51 

69.13335 

66 

69.29471 

52 

69 . 14529 

67 

69.30361 

53 

69.15711 

68 

69.31219 

54 

69 . 16882 

69 

69.32046 

55 

70 

We  shall  now  apply  the  above  projection  to  the  construction  of  a 
map  of  the  Dominion,  by  assuming 

<//  =  42°  </>"  =  70° 

and  the  extreme  longitudes  to  be  55°  and  141°  west.  The  latitudes  of 
the  standard  parallels  are  also  assumed  to  be: 

0i  =  45°,  02  =  60°. 

By  eq.  (10),  and  the  above  table,  we  find 

fi  =  3548.3829m 
r2  =  2511. 2187m 

Also,  by  eq.  (13) 

n  =  0.791123. 

Also,  by  eq.  (14) — * 

<r'=  1.00722,  for  lat.  42° 

<r'=  1.05811,  for  lat.  70°. 

These  are  also  the  area  scales  at  the  extreme  latitudes. 

The  latitude  of  the  parallel,  between  the  two  standard  parallels,  at 
which  the  scale  ratio  in  longitude  is  a  maximum  is,  eq.  (18): 

0  =  52°  56'  05" 

and  for  this  latitude 

o'  =  0.9914455 
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For  this  latitude  we  also  find — eq.  (17): 

a'-a=  =  -14'  46". 04 

At  latitude  63°  we  find 

<r'=  1.010086 

so  that  up  to  latitude  63°  the  scale  error  will  barely  exceed  one  per  cent. 

By  assuming  the  latitudes  of  the  standard  parallels  differently  a 
different  distribution  of  scale  error  may  be  obtained. 

(2)  The  maximum  scale  error  between  the  standard  parallels  is  to 
be  equal,  with  opposite  sign,  to  that  of  the  extreme  parallels. 

Let  z  denote  the  length  VP,  P  being  the  projection  of  the  pole,  and 
m'  and  m"  the  lengths  of  the  meridian  arcs  between  the  extreme  parallels 
and  the  pole. 

Then,  the  scale  ratios  of  these  parallels  being  equal,  we  have 

(z-\-m')0  ( z-\-m")6 


whence 


P'  CO 

m'P " 


P" 


CO 


...ii  T>r 

m  P 


z  — 


(20) 


pr _ p " 

which  determines  z. 

Again,  to  find  the  latitude  corresponding  to  the  maximum  scale  error 
between  the  standard  parallels,  we  have,  for  a  parallel  near  the  centre 
of  the  map 

nizpm), 


l-a'=l 


or 


<T  = 


P 

n(z-\-m) 


(21) 


We  now  differentiate  this  with  reference  to  <f>  and  equate  to  zero,  thus 

dm  n  dP 

P  —  —  (z+ra)  — 

da'  deb  deb 

—  =n - =  0 


d(p 


P2 

dP 


Then  as 
this  becomes 


-rj  dm  ,  % 

P  —  =  (z  +  m) 
d<p  d<f> 

dm  ,  dP 

— -  —  —  p  and  — -  =  —  p  sin  $ 

d4>  d(f> 


Pp=  (z+w)p  sin  </>, 

or  N  cos  <f>  =  (zpm)  sin  </>, 

or  Ao  cot  0O  —  m0  =  z, 

which  gives  the  latitude  corresponding  to  the  maximum  scale  error 
This  equation  must  be  solved  by  trials. 


(22) 
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Again,  the  scale  ratio  being  the  same,  with  opposite  sign,  on  this 
parallel  and  either  extreme  parallel,  we  have 

n(z-\-m0)  _  n(z-\-m') 


1- 


P  r 


P' 


-1 


whence 


1 

n 


z+m0 


z-\-mf 
2  P' 


which  gives  n. 

Again,  for  the  standard  parallels 

n(z-\-nii)  _  n(z-\-ni2 ) 
P2 


or  P  \  —  nmx  —  nz,  (24) 

and  P2  —  wm2  =  nz; 

by  which  0i  and  02  may  be  determined  by  trials. 

This  modification  of  (1)  may  now  be  applied  to  the  construction  of 
a  map  of  Canada  by  assuming  again  the  extreme  latitudes  to  be 

0'  =  42°  (/>"  =  70° 


We  must  first  find  m'  and  m"  By  means  of  eq.  (8)  we  find — 

m"=  1387^5858 

and  by  adding  to  this  the  lengths  of  degrees  of  latitude  contained  in 
the  table,  p.  1 1  we  find- — -  m'  —  3324w6545. 

Then  by  table  in  D.L.S.  Manual — - 

log  P'  =  3.4697861, 
log  P"  =  3. 1334070. 

Then  eq.  (20)  gives 

2  =  268.5985. 

Then  by  eq.  (22) — - 

(/>o  =  57°  49'  20". 17. 

The  value  of  n  then  follows  by  (23) : 

n  =  0.8334627.  .  .  ; 
for  which  the  scale  ratio  is  found  by  (21): 

a'  =  0.9847151 
=  1-0.0152849 

As  a  check  the  scale  ratio  is  found  for  either  extreme  parallel  to  be 

a' =1.0152848 

The  maximum  error  in  a  direction  for  latitude  </>0  is  found  by  (17)  to  be 

a'-a=  -26'  28". 5; 

and  this  is  the  maximum  for  the  whole  map. 

The  latitudes  of  the  standard  parallels  are,  by  (24) 

</>!  =  46°  53'  39".  15 
<£2  =  660  53'  34". 81. 

(3)  Albers’  Equivalent,  or  Equal-Area,  Projection. 
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A  brief  development  of  the  theory  of  this  projection  may  not  be  out 
of  place. 

The  scale  ratios  in  latitude  and  longitude  may  be  thus  expressed 

nr 


a  — 


dr  ,  rdO 
-  <j  = 


pdcp  Pda,  N  cos  0 

Also  the  condition  for  equivalence  of  areas  is: 

a  a'  =  1 

dr  nr 


(25) 


or 


=  1 


pd cp  N  cos  0 

or  nrdr  =  —  pN  cos  (pd(p 

By  (1)  and  (2)  this  becomes 

,  a2(  1-e2) 

nrdr  =  — - cos  <pd(p 


(26) 


or 


rdr  =  — 


(1  —  e2  sin2  0)2 
a2(l—e2)  cos  <pd<p 


n 


(1  —  e2  sin2  (py 


Therefore 


j/2 


a2(l  —  e2) 


cos  <pd<p 


(l—e2  sin2  4>y 


2  n 

If  now  r0  denote  the  value  of  r  when  (p  =  0  this  becomes 

0  cos  (pdcp 


,  ,  2a2  (1  —e2) 

r~  -  r0“  =  -  - 

n 


o  (1  —e2  sin2  (p)2 


If  then  the  earth  were  a  sphere  having  the  radius  R ,  and  13  denote  the 
value  of  the  latitude  on  this  supposition,  the  right-hand  member  of 
the  last  expression  would  become: 

2R2  'I3  RJR  2R2  •  « 

-  cos  (3d(3=  — -  sin  j3. 

n  Jo  n 

We  may  therefore  define  (3  by  equating  the  two  expressions,  which  gives 

7")o  •  o  2/i  <>\  0  cos  <pd(p 

R~  sin  /3  =  U“(1 — e2)  - — 

Jo  (1  —  e2  sin2  <p)2 

Therefore,  expanding  and  integrating  each  term  separately,  we  have 

’0 


R2  sin  jS  =  a2(l  —  e2 


(l  +  2c2  sin2  0-fi- 3c4  sin4  0 
+  4e6  sin6  0  +  )cos  (pdcp 


9  p2 

=  a2(  \—e2)  (sin  0+ —  sin3  0-f-  —  sin5^> 

3  5 

4c6 

+  —  sin70  +  ). 

7 

•  •  7T  7T 

R  may  be  determined  by  assuming  that  when  0  =  -  ,  f3  also  =  -,  which 

2  2 
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gives 


R2  =  a2(  1-e2)  (  l+?l2  +  —  +  —  + 

3  5  7 


(27) 


R  is  the  radius  of  the  sphere  having  the  same  surface  area  as  that 
of  the  spheroid.  For  the  Clarke  spheroid  of  1866 — 

log  R(in  metres)  =6.80420742 
log  i?(in  miles)  =3.59756243 
/3  is  then  given  by  the  expression: 


sin  j8  =  sin  0 


-  .  2e2  .  3e4  .  4e6  . 

1+  — -  sin-  0  +  —  sin4  0-f  —  sin6  0  + 

3  5  7 

i  +  2A!  +  3_£.4  +  ii6  + 

3  5  7 


=  sin  0 


\ 


3 

64 


cos 20  — (  —  cos20  +  —  sin2  0  cos2  <j>  )e 4 
\45  5 


6  4 

cos2  0+  —  sin2  0  cos2  0  +  —  sin4  0  cos2  0  )e6 
,945  35  7  ' 


+  etc. 


(28) 


is  termed  the  “authalic”,  or  equal-area,  latitude.  The  difference 
0  —  (3  may  be  computed  by  the  expression: 

0-/3  =  467". 0129  sin  20-O".4494  sin  4 0 

+  0".0005  sin  6 0 — -  (29) 

A  table  is  given  in  Special  Publication  No.  67  of  the  U.S.  Coast  and 
Geodetic  Survey  which  contains  values  of  0  —  (3  for  30'  intervals  of 
latitude  from  0°  to  90°, 

The  expression  for  r  may  then  be  written: 


r-rQ- 


2  2  R*  . 


sin  (3 


n 


(30) 


The  constants  n  and  ro  must  now  be  determined.  Introducing  the 
condition  that  the  scale  shall  be  true  along  the  standard  parallels,  we 
have 


N i  cos  <f>i  7V2  cos  </)2 
r  i  = -  r2  =  - 


(31) 


n  n 

and  substituting  these  values  of  r\  and  r2  in  the  expression  for  r  we  have 


2  R2  . 


rl - sin  /3i  = 


N\  cos2  0i  P\ 


n 

2  R2 


n 2 


ro~ 


sin  jS2  = 


N\  COS2  02  Pj 
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Then  subtracting  we  find 

2  R2 


sin  /32  —  sin  /3i  )  =  — — — 


n 


whence 


n  = 


P2  —  P2 
1  1  r  2 


2P2  (sin  jS2  —  sin  (fi) 

Again,  from  the  general  expression  for  r  we  have 


(32) 


2  2  2  R2  . 

rx  =  r0 - —sin  ft 

n 

whence 

2  2P“  (  ‘  Q 

r-  —  rx—  -  (sin  ft  — sin  p) 

n 

or 

OR2 

P  —  r]- 1 - (sin  /Si  —  sin  ft. 

n 

(33) 

Similarly 

2R2 

P  —  r\-\ - (sin  ft  —  sin  ft 

n 

(34) 

For  convenience  of  computation  the  expressions  for  r  may  be  written 

(35) 


-  =(  '-L  +  -  sin  13,)-  2  sin  0 
R 2  \R2  n  )  n 


T2  _[y\ 


r 2  (r\  2  .  „  \  2 


-  =  (  -  +  -  sin  & 
R 2  \P2  n 

The  scale  ratio  for  any  parallel  is 


sin  j3 


n 


(36) 


(7  = 


nr 

P 


(37) 


and  that  for  a  meridian  at  any  point 

1 


P 


(38) 


a  nr 

To  find  the  latitude  of  the  parallel,  between  the  standard  parallels, 
for  which  a  is  a  maximum,  we  have 

d'  /  Tv  T  \  AT  dr 

r  —  (A  cos  <f>)  —  A  cos  </>  — 

dcr  1  dcf)  d(f> 

d<f>  n  r2 

so  that,  equating  to  zero,  the  condition  for  a  maximum  is 

r  —  (A  cos  <j>)  =  N  cos  (f>  —  . 
dcf>  d4> 

TA  d  /AT  \  -T  . 

But  —  (A  cos  </>)  =  cos  <£  —  —  A  sin 
d<f>  dcf) 


Engineering  Research  Bulletin 


17 


a  -A 

=  —  cos  <j>  —  (1  —  e2  sin2  <f>)  j(  —  2e2  sin  0  cos  0) 
—  A7"  sin  0 

ac2  sin  0  cos2  0  —  a  sin  0  (1—e2  sin  20) 

(1  —  e2  sin2  0)"^ 
a(l  —  e2)  sin  0 

=  —  - -  =  -  p  sin  0 

(1  —  e2  sin2  0)  ^ 

Also,  by  (26)  *  -  -  C0S  *  , 

d(j)  nr 

Therefore  the  condition  for  a  maximum  becomes 

N  cos  4>pN  cos  0 


fp  sin  0  = 


nr 

pN 2  cos2  0 


wr 

79 


or 


N2 

tan  0  sec  0  = - 

nr1 


(39) 


By  eq.  (16)  we  have  for  the  maximum  error  in  direction  at  a  point 

1 


a/  —  a  = 


1  —  (T2 

(1T<t2)  sin  1" 


(40) 


We  now  apply  this  projection  to  the  construction  of  a  map  of  the 
Dominion,  by  assuming  as  before 

0'  =  42°  0"  =  70° 

We  shall  also  take  the  standard  parallels  in  the  latitudes 

4>i  =  45°  02  =  65° 

By  eq.  (32)  we  find 

n  =  0.8067798  .  .  . 

Also  by  eq’s,  (31) 

n  =  3479.5229 
r2  =  2081.8885 

The  latitude  for  which  <j  is  a  maximum  is  found  by  a  series  of  approxi¬ 
mations,  using  eq.  (39),  to  be: 


56°  17'  05"  .04; 

for  which  we  find 

a  =1.01539135  .  .  . 
Also  <r'  =  0.9848418  .  .  . 

=  1-0.0151582.  .  .  . 
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We  also  find — eq.  (40) — -for  the  same  latitude: 

a'  — a=  -  52'  30". 28. 

The  following  table  gives  the  value  of  r  for  each  degree  of  the  map 
computed  by  eq.  (36),  from  which  are  found  the  length  of  each  degree 
of  the  developed  meridian.  The  errors  of  the  latter  are  found  by  com¬ 
parison  with  the  lengths  given  in  the  table,  p.  11. 


Lat. 

r 

Length 

Error 

of  1° 

O 

m 

m 

m 

or 

Yo 

42 

3685.8212 

68.5672 

-0.4563 

-0.66 

43 

3617.2540 

68 . 7695 

-0.2662 

-0.39 

44 

3548.4845 

68.9616 

-0.0364 

-0. 13 

45 

3479.5229 

69.1430 

0.0827 

0.12 

46 

3410.3799 

69.3137 

0.2411 

0.35 

47 

3341  0362 

69.4725 

0.3877 

0.53 

48 

3271.5937 

69.6193 

0.5223 

0.76 

49 

3201.9744 

69.7530 

0.6438 

0.93 

50 

3132.2214 

69 . 8734 

0.7521 

1.09 

51 

3032.3480 

69.9789 

0.8453 

1  22 

52 

2992.3691 

70.0393 

0.9243 

1.34 

53 

2922 . 2995 

70. 1438 

0.9857 

1.43 

54 

2852.1557 

70.2000 

1.0312 

1.49 

55 

2781  9557 

70.2380 

1 . 0576 

1.53 

56 

2711.7177 

70.2552 

1 . 0334 

1.54 

57 

2641.4625 

70.2513 

1 . 0482 

1.51 

58 

2571  2112 

70.2226 

1 . 0385 

1.46 

59 

2500.9886 

70. 1684 

0.9434 

1.36 

60 

2430.8202 

70.0S31 

0.8504 

1.23 

61 

2360.7341 

69.9725 

0.7263 

1.05 

62 

2290.7616 

69.8244 

0.5680 

0.82 

63 

2220.9372 

69.6384 

0.3720 

0.54 

64 

2151.2988 

69.4103 

0. 1342 

0.19 

65 

2081.8885 

69.1349 

-0. 1506 

-0.22 

66 

2012.7536 

68 . 8062 

-0.4885 

-0.70 

67 

1943.9474 

68.4185 

-0.8851 

-1.28 

68 

1875.5289 

67 . 9625 

- 1 . 3497 

-1.95 

69 

1807.5664 

67.4307 

- 1 . 8898 

-2.73 

70 

1740. 1357 

(4)  The  Lambert^Conformal  Projection. 

A  brief  development  of  the  theory  of  this  projection  is  also  given  here. 
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F/g.  2 


In  Fig.  2  A  BCD  represents  an  elementary  quadrilateral  on  the 
surface  of  the  earth,  bounded  by  meridians  and  parallels  of  latitude, 
and  abed  is  its  representation  on  the  map. 

If  u  denotes  the  co-latitude  of  AB  we  have 

a(  1  —  e2)  du 


BD  =  pdu  = 


AB  —  N  sin  udoo  = 


(1  —  e2  cos2  u) f 

a  sin  udoi 


( 1  —  e2  cos2  u) 5 

Then,  as  the  representation  is  conformal 

bd  _  BD 
ab  AB 

dr  a(l  —  e2)du  (1— e2cos 2  u)\ 


or 


nrdoi 


(1  —  e2  cos2  w)-f  a  sin  ndco 
1  dr  _  (1—  e2)du 

n  r  sin  u(l  —  e2  cos2  u) 


1 


sin  u 
du 


1 


2  1  —  COS2  U 

1—e2  cos2  ut 
e 2  sin  udu 


\du 


sin  u  1—  e2  cos2  u 

Then,  resolving  into  partial  fractions,  this  becomes 

1  dr  du  e2  sin  udu  e2  sin  udu 


2(1  +  0  cos  u)  2(1  —  e  cos  u) 


(41) 


n  r 


sin  u 
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Then  integrating,  we  have 

u  e 

—  log  r  =  log  tan  —  +  —  log  ( 1  +  e  cos  u) 
n  2  2 

6  1 

—  -  log  (l  —  e  cos  u)-\-  —  log  k 

2  n 

the  last  term  being  the  constant  of  integration.  This  may  be  written 

.  r  ,  i  u  ( l-\-e  cos 

log—  log |  tan  -  1  -  12 


or 


r=k[  tan 


u 


7i+ 


l  —  e  cos  u, 

ne 


6  COS  U 


2/  \1- 


e  cos  Uj 


If  now  an  angle  z  be  assumed,  such  that 

z  u  (  l-\-e  cos  u\~ 

tan-  =tan  —l  -  12 

2  2  V 1  —  e  cos  u) 

and  also  an  angle  q,  such  that 

cos  q  =  e  cos  u 

we  shall  have 


so  that 


1  +  e  cos  u 


e  cos  u 


cot  -  )  = 


1  +  cos  q 
,1  —  cos  q, 

1  +  c  cos  u\~ 
l  —  e  cos  u) 


=  cot  |  q; 


(42) 


and 


Z  U 

tan  —  =  tan  —  (  cot 


n 


r = k ( tan 


(43) 

(44) 


A  close  approximation  to  the  angle  z  may  be  found  by  computing 
the  geocentric  latitude  <j>r  by  the  equation: 

tan  <//  =  (1  —  e2)tan  4>;  (45) 

then  approximately 

2  =  90°  —  (//  (46) 

The  value  of  n  may  now  be  found  by  assuming  the  scale  to  be  true 
along  the  two  selected  standard  parallels.  Then  making  use  of  the  fact 
that  the  lengths  of  two  arcs  of  these  parallels,  each  having  the  longitude 
co,  are  proportional  to  their  representations,  we  have 

WCOfi  _  Ni  COS  (f)  1  .  CO 
WC0f2  N2  COS  c/>2  .  CO 

r\  N\  cos 

or  —  =  - 

r-2  N2  cos  02 
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or 


tan  \zy 


N i  cos  Pi 


dan  \zj  N%  cos  02  P2 
Then  taking  log’s,  and  solving  for  n,  we  have 

log  Pi -log  P2 


n  = 


1  2i  .  S2 

log  tan  —  —  log  tan  - 
2  2 


(47) 


An  expression  for  k  may  now  be  found  by  equating  the  length  of  an 
arc  of  a  standard  parallel  and  that  of  its  representation,  thus 

nr  ico  =  Pico 


or 


or 


nk[  tan  —  )  =Pi 


k  = 


Pi 


P‘1 


Si 


w(  tan 

V  2 


w 


wl  tan  — 


(48) 


The  collected  equations  then  are: 

cos  q  =  e  cos  u  —  e  sin  0 


2  u  q 

tan  -  =  tan  -  I  cot  - 

2  2  V  2 


=  tan  (45°  — 10)  (  cot  — 


n  = 


log  Pi  — log  P2 

log  tan  —  —log  tan  — 

2  2 


k  — 


Pi 


P2 


n  (  tan  — 
2 


r  =  k[  tan  - 


/1 


w^tan  ^ 


n 


(49) 


Eq’s.  (19)  may  be  used  in  plotting  a  parallel  ofjatitude  by  coordi¬ 
nates. 

The  scale  ratio  of  any  parallel  is  evidently: 

,  nwr  nr 

= - —  7 

Pco  P 

see  eq.  (14).  Also,  as  the  projection  is  conformal 

/ 

U  =  <7 

The  area  scale  is 

'=  ij™  y. 

VP  /  ’ 

Eq,  (16)  also  gives  a'  —  a  =  0; 

so  that  directions  are  correctly  represented  at  a  point. 


(50) 

(51) 

(52) 
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To  find  the  latitude  for  which  the  scale  error,  between  the  standard 
parallels,  is  a  maximum,  we  have,  as  on  p.  10  the  necessary  condition: 


Then,  by  (41) 


dr 

d(f) 


=  nr 


-.dr  dP 

P  =  r— 
d4>  d<f> 


e 2  cos  0 


1 


1  —  e-  sin2  0  cos  <f>, 


Also  on  p.  10  it  was  shown  that 

dP 


d(f> 


=  —  p  sin  0. 


Therefore  the  above  condition  becomes 

cos  *  1 

a  cos  0 


1  —  e~  sin2  <(>  cos  0, 


=  —  p  sin  0 


or 


or 


n 


e 2  cos2  0  —  1  +  e2  sin2  4> 

-  =  —  p  sin  4> 


( 1  —  e2  sin2  4>) 2  cos  0  ( 1  —  e2  sin2  0) 

o(l  —  e2) 


—  n 


=  —  p  sin  0 


(53) 


(1  —  e2sin20)f 

or  np  =  p  sin  0 

or  sin  0  =  w, 

the  required  condition. 

Applying  this  projection  to  the  same  data  as  in  the  previous  case,  viz. : 

0'  =  42°,  0"  =  70° 

0i  =  45°,  02  =  65° 
we  obtain  the  following  results: 

Eq.  (47)—  n  =  0.8234698  .  .  . 

Eq.  (31)—  n  =  3409.0005 

r2  =  2039.6931. 

Also  the  latitude  for  which  the  scale  error  is  a  maximum  is — -eq.  (53) : 

0  =  sin-1w  =  55°  26'  01".16; 
and  for  this  latitude  we  find — -eq.  (49)  : 

=  =  0.984733  .  .  . 

=  1-0.015  .  .  . 

Also  the  area  scale — -eq.  (51): 

02  =  0.969699 
=  1-0.030.  .  .  . 

Directions  are  correctly  represented — -eq.  (52). 

The  properties  of  projections  (3)  and  (4)  are  evidently  mutually 
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exclusive,  as  it  is  impossible  for  any  projection  to  embody,  for  all  parts 
of  the  map,  the  two  conditions: 

( Tv'  =  1,  and  a  =  ar , 

unless,  at  the  same  time,  it  were  possible  that — - 

<7=1,  and  a'  =  1. 

This  could  only  be  realized,  however,  if  the  surface  of  a  sphere  or  a 
spheroid  were  developable. 

The  following  table  gives,  for  Lambert’s  projection,  the  value  of  r 
for  each  degree  of  latitude  between  42°  and  70°,  from  which  the  lengths 
of  the  individual  degrees  of  latitude  follow.  The  errors  of  the  latter 
are  found  by  comparison  with  the  lengths  given  in  the  table,  p.  11. 


Lat. 

r 

Length  of  1° 

Error 

O 

m 

m 

m 

% 

42 

3617.0780 

69.5744 

0.5509 

0.80 

43 

3547 . 5036 

69.3532 

0.3175 

0.46 

44 

3478.1504 

69.1499 

0.1019 

0.15 

45 

3409.0005 

68 . 9644 

-0.0959 

-0.14 

46 

3340.0361 

68.7972 

-0.2754 

-0.40 

47 

3271.2389 

68.6478 

-0.4370 

-0.63 

48 

3202.5911 

68.5166 

-0.5804 

-0.84 

49 

3134.0745 

68.4040 

-0.7052 

-1.02 

50 

3065 . 6705 

68.3096 

-0.8117 

-1.17 

51 

2997 . 3609 

68.2341 

-0.8992 

-1.30 

52 

2929 . 1268 

68. 1775 

-0.9678 

-1.40 

53 

2860.9493 

68.1402 

-1.0169 

-1.47 

54 

2792.8091 

68.1227 

-1.0461 

—  1  51 

55 

2724 . 6864 

68.1251 

-1.0553 

-1.53 

56 

2656.5613 

68.1482 

- 1 . 0436 

-1.51 

57 

2588.4131 

68.1924 

-1.0107 

-1.46 

58 

2520.2207 

68.2586 

-0.9555 

-1.38 

59 

2451.9621 

68.3474 

-0.8776 

-1.27 

60 

2383.6147 

68 . 4593 

-0.7764 

—  1  12 

61 

2315.1554 

68.5956 

-0.6506 

-0.94 

62 

2246.5598 

68.7576 

-0.4988 

-0.72 

63 

2177.8022 

68 . 9463 

-0.3201 

-0.46 

64 

2108.8559 

69.1628 

-0.1133 

-0.16 

65 

2039.6931 

69 . 4090 

0.1235 

0.18 

66 

1970.2841 

69.6870 

0.3923 

0.57 

67 

1900.5971 

69 . 9986 

0.6950 

1.00 

68 

1830.5985 

70.3463 

1.0341 

1.49 

69 

1760.2522 

70.7327 

1 . 4122 

2.04 

70 

1689.5195 

The  Polyconic  Projection 

In  this  projection  each  parallel  of  latitude  is  assumed  to  be  the  base 
of  a  right  circular  cone,  tangent  to  the  surface  of  the  earth  along  the 
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parallel,  and  whose  vertex  is  consequently  on  the  axis  of  the  earth 
produced.  When  the  cone  is  developed  the  parallel  becomes  a  circular 
arc  whose  radius — • 

r  =  N  cot  0,  (54) 

The  central  meridian  of  the  map  is  represented  by  a  straight  line,  on 
which  the  true  lengths  of  the  degrees  of  latitude  are  laid  off  to  scale. 
Each  parallel  is  then  described  about  its  own  centre,  and  the  degrees 
of  longitude  laid  off  on  it  correctly  to  scale.  Through  the  points  of 
division  thus  found  the  remaining  meridians  are  drawn, 

As  the  scale  is  true  along  any  parallel  we  have  the  relation: 

rd  =  Eco ; 

or  6N  cot  0  =  coiV  cos  0; 

so  that  6  =  oo  sin  </>;  (55) 

0 

and  n—  —  =  sin  0.  (56) 

CO 

n  has  therefore  a  different  value  for  each  parallel. 

The  parallels  may  be  drawn,  on  maps  of  large  scale,  by  means  of 
coordinates  computed  by  (19). 

The  scale  of  the  map  is  thus  correct  along  the  central  meridian, 
and  along  the  parallels,  but  the  scale  in  latitude  increases  considerably 
with  increase  of  distance  from  the  central  meridian.  The  meridians, 
in  general,  are  represented  by  curved  lines,  and  they  do  not  intersect 
the  parallels  at  right  angles. 

The  increase  of  scale  in  latitude  is  given  approximately  by  the 
expression  (See  U.S.C.  and  G.  Survey  Spec.  Pub.  No,  57) : 

E  =  0.0l(G^iiY  (57) 

in  which  co  is  expressed  in  degrees.  By  means  of  this  equation  we  find 
the  scale  error  at  some  points  as  follows: 

0  =  49°,  Z  =  125°:  £  =  0.0410  or  4.10%. 

0  =  5 5°  26',  Z  =  130°:  £  =  0.0442  or  4.42%. 

0  =  70°,  Z  =141°:  £  =  0.02997  or  2.997%. 

Let  us  now  find,  for  latitude  49°,  the  range  of  longitude  for  which  the 
scale  error  in  latitude  amounts  to  0.01.  Placing  £  =  0.01  in  (57)  and 
solving  for  co  we  find 

co  =  8.1  sec  0 
=  12°. 346  .  .  . 

for  latitude  49°.  This  projection  should  therefore  not  be  used  for  maps 
covering  a  wider  range  in  longitude  than  about  25°.  It  is,  however, 
within  these  limits,  a  very  useful  projection,  on  account  of  the  ease 
with  which  it  can  be  constructed;  extensive  tables  being  available  which 
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are  applicable  to  any  part  of  the  world,  and  which  greatly  facilitate 
the  work  of  the  cartographer. 

Bonne’s  Projection 

This  is  a  modification  of  the  simple  conical  projection.  The  scale 
is  true  along  a  central  parallel,  and  the  degrees  of  latitude  are  laid  off 
correctly  along  the  central  meridian  which  is  represented  by  a  straight 
line.  The  radius  of  the  central  parallel: 

r0  =  No  cot  <t>0,  (58) 

<t> o  being  its  latitude,  and  the  remaining  parallels  are  circular  arcs  des¬ 
cribed  about  the  same  centre,  so  that  the  radius  of  any  parallel: 

r  =  r0±m,  (59) 

m  being  its  distance,  measured  along  the  central  meridian,  from  the 
central  parallel.  The  degrees  of  longitude  are  then  laid  off  truly  along 
each  parallel  and  the  points  of  division  joined,  as  in  the  polyconic  pro¬ 
jection.  We  have  then,  for  any  parallel,  the  relation 

r6  —  N  cos  (/> .  co 
N 

whence  6=  —  co  cos  <f>  (60) 

r 

The  points  of  intersection  of  the  meridians  and  parallels  may  be 
plotted  by  means  of  coordinates,  computed  by  the  equations: 

x  =  r  sin  6 

y  =  ro  —  r  cos  6  (61) 

the  origin  being  the  point  of  intersection  of  the  central  meridian  and 
the  central  parallel. 

A  point  in  favour  of  this  projection  is  that  it  is  equal-area;  the  main 
objection  to  it  is  on  account  of  the  obliquity  of  the  intersections  of  the 
meridians  and  parallels  at  points  distant  from  the  central  meridian. 
The  divergence  from  a  right  angle  at  any  such  intersection  may  be 
found  by  the  expression 

tan  \p  =  {0  —  co  sin  0),  (62) 

180 

d  and  co  being  expressed  in  degrees.  For  latitude  70°  and  longitude  141° 
(co  =  41°)  this  gives  \p  =  5°  54'  20". 

Also  the  scale  error  in  latitude  at  this  point  is 

0.005  or  0.5  per  cent. 

At  the  point  whose  geographical  coordinates  are 

</>  =  49°,  L  =  125° 
the  scale  error  in  latitude  is  zero,  and 

^  =  1°  43'09".6. 
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The  scale  errors  in  latitude  thus  appear  to  be  small.  This  advantage, 
however,  is  offset  by  the  fact  that  the  obliquity  of  the  intersections  of 
the  meridians  and  parallels  causes  the  two  diagonals  of  a  quadrilateral 
formed  by  two  meridians  and  two  parallels,  which  should  have  the 
same  length,  to  have  widely  different  lengths.  This  fault  does  not 
exist  in  the  four  modifications  of  the  conical  projection  with  two  standard 
parallels,  considered  above,  in  which  the  intersections  are  all  right- 
angled. 


